Tracing magnetism and pairing in FeTe-based systems 
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In order to examine the interplay between magnetism and superconductivity, we monitor the non- 
superconducting chalcogenide FeTe and follow its transitions under insertion of oxygen, doping with 
Se and vacancies of Fe using spin-polarized band structure methods (LSDA with GGA) starting 
from the collinear and bicollinear magnetic arrangements. We use a supercell of FesTes as our 
starting point so that it can capture local changes in magnetic moments. The calculated values 
of magnetic moments agree well with available experimental data while oxygen insertions lead to 
significant changes in the bicollinear or collinear magnetic moments. The total energies of these 
systems indicate that the collinear-derived structure is the more favorable one prior to a possible 
superconducting transition. Using a 8-site Betts-cluster-based lattice and the Hubbard model, we 
show why this structure favors electron or hole pairing and provides clues to a common understanding 
of charge and spin pairing in the cuprates, pnictides and chalcogenides. 

PACS numbers: 65.80. +n, 73.22.-f, 71.27.+a, 71.30.+h 



High temperature superconductors discovered in the 
1980s consisted of copper-oxide-based, layered materi- 
als. In such cuprates, it is believed that doping, away 
from half filling, of a Mott-insulator leads to supercon- 
ductivity although there is still no general agreement on 
a specific mechanism. Recently discovered superconduc- 
tivity in Fe-based compounds, which are either pnictides 
or chalcogenides [U |2] , has opened up an extremely rich 
and active area of basic research. Although the transition 
temperatures of the Fe-chalcogenide superconductors are 
among the lowest of the recently discovered compounds, 
they possess rather simple layered structures and fas- 
cinating antiferromagnetic or spin density wave (SDW) 
states. Superconductivity in the Fe-chalcogenides, hav- 
ing the so-called (11) structure, was first reported in 2008 
(Refs. 013] ) in Fei+^Se and Fe(Sei_^Te^)o.82. This dis- 
covery led to a substantial increase in research efforts 
focused on simple, layered Fe compounds, containing 
chalcogenides such as S, Se, and Te. FeTe, which is a 
metallic antiferromagnet, has a tetragonal structure and 
has shown properties uniquely different from some of the 
other pnictides and chalcogenides. For example, Fei+^Te 
is not superconducting; instead, it shows magnetic and 
structural transitions at 65 K [4 . In addition, it is said 
to have the so-called "double stripe" antiferromagnetic 
order instead of the "single stripe" with ordering vector 
(0.5, 0.5) [5 . Nevertheless, FeTe doped with Se was found 
to be superconducting (Refs. [HE]). Oxygen insertions 
in FeTe (films) have also given rise to superconductivity 
(Refs. [TjISl). There are claims that the Fe-chalcogenides 
do not exhibit the same nesting feature found in the pnic- 
tides and that they possess a different local SDW mag- 
netic order which survives even in the highest Tc sam- 
ples [9]. 

The present work is focused on understanding such 
changes while attempting to find common behavior, es- 



pecially en route to superconductivity, in these chalco- 
genide (and possibly pnictide and cuprate) systems. We 
note that there have been numerous band structure cal- 
culations carried out recently on related systems focusing 
on nesting features of the relevant Fermi surface, suscep- 
tibilities among other things (see Ref. [10^ and references 
therein); however, to our knowledge, these studies have 
not focused on systematic changes in the resulting mo- 
ments and/or total energy as functions of various modu- 
lations en route to superconductivity. 



Band structure 

Our band structure calculations are based on the 
(spin-polarized) density functional theory (DFT) [TTJ 
12 , combined with the generalized gradient approxima- 
tion (GGA) [13], as implemented in the VASP pack- 
age [T4HT6] . The interaction between the electrons and 
atomic cores is described by projector augmented wave 
(PAW) pseudopotentials. The wave functions are ex- 
panded using plane waves with a cut of energy of 400 
eV. Brillouin Zone integrations are carried out with a 
4x4x4 Monkhorst-Pack [17] grid of k-points while a 
denser 8x8x8 k-point mesh was used for the density 
of states (DOS) calculations. The lowest energy struc- 
ture was determined using Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) [18 based algorithms. For the results 
reported in this work, all the atoms are allowed to re- 
lax based on Hellman-Feynman forces and optimization 
is carried out until typical forces on the atoms are around 
0.01 eV/A or less. 

FeTe has a tetragonal structure with space group 
P4/nmm, where a square lattice of Fe is tetrahedrally co- 
ordinated with Te ions. Experimental lattice parameters 
of the pure FeTe are a = 6 = 3.8248 A and c = 6.2910 A; 
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FIG. 1: Projected FeTe density of states due to various mod- 
ulations and the supercell studied here; Fe atoms are located 
on the plane shown and possible oxygen interstitial sites are 
identified as A, B, C and D; Te atoms are gray spheres in 
between the planes of Fe. 

these were used as the starting lattice pararaeters for the 
2x2x1 supercell. This 16-atom tetragonal supercell, 
consists of 8 Fe and 8 Te atoms (see Fig. [ijb)) with Fe 
atoms located in a plane with Te atoms sitting above 
(and below) the plane. The reason for choosing a large 
supercell is to examine possible nontrivial SDW-type be- 
havior; i.e., to allow moments to (vary and) align among 
themselves as necessitated by the variational principle 
in total energy. A smaller cell will not have such free- 
dom. Following the (non-spin-polarized) work of Ref. , 
the interstitial sites for oxygen insertion (A, B, C, D in 
Fig. [TJb)) were selected by comparing the total energies 
of several possible interstitial sites containing oxygen in 
the unit cell. Although our spin-polarized work shows 
some differences (related to sites C and D) with the above 
unpolarized work, sites A and B (for oxygen interstitials) 
appear to be the most energetically favorable, as deter- 



mined there. 

The other modulations studied in this work are; (a) 
Fe interstitials, (b) Fe vacancies and (c) Se dopants. It 
is experimentally known that, except in case (a), these 
changes to pure FeTe, carried out with sufficient care, are 
likely to induce superconductivity. We note that without 
the GGA, the calculated moments turn out be rather 
small and (quantitative) conclusions drawn from pure 
LDA (or LSDA) work alone may be questionable. 

Magnetic Moments 

The magnetic moments in pure FeTe and its various 
modulations were calculated starting from two different 
(AFM) magnetic structures, collinear and bicollinear (see 
Fig. [5]). Note that in the collinear structure, k = (0.5, 0.5) 
(with respect to a 2-Fe unit cell) antiferromagnetic or- 
der is present while the bicollinear structure has (0.5, 
0.0) order. The pure compound FeTe, in both mag- 
netic structures, carry Fe moments which are around 2.0 
jiiB, in agreement with recent experiments. This fact 
alone points to the relevance and importance of multi- 
orbital effects at the Fermi surface since otherwise the 
moments cannot be this high. Our density of states plots 
(Fig. [ija)) identify these as originating from d-orbitals, 
which is not a surprise. These moments are found to be 
severely affected by the insertion of two oxygen atoms 
into interstitial sites A and B (see Table |l] and [ll| ; Four 
of the eight Fe atoms show moments that are around 1.2 
/i^ or less while the total charge in each atom under- 
goes minimal (less that one percent) changes. In den- 
sity functional calculations based on various approxima- 
tions such as the LDA or GGA, such changes are not un- 
usual. Hence, one cannot make specific statements about 
changes in the valence (such as Fe^+ Fe^+) using these 
results alone. However, when the total valence charges 
of Te and O are compared (at the same radii), there is 
a clear difference. Oxygen has more charge while being 
more electronegative and hence, its insertion is likely to 
create holes in the Fe atoms. 

Insertion of Fe into A and B interstitial sites of FeTe 
does not appear to change the moments as drastically. In 
fact, while the moments of the Fe atoms in FeTe are al- 
most unchanged, showing long-ranged antiferromagnetic 
order, the inserted Fe atom carries a moment of about 
2.46 /iB, close to the value reported by some experi- 
ments [19 . The variationally determined total energies 
of these two, collinear and bicollinear, structures show 
a trend that supports experimental findings for FeTe 
and Se-doped FeTe. With the pure system exhibiting 
collinear (static) antiferromagnetic (AFM) ground state. 
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FIG. 2: Damped antiferromagnetic or SDW arrangements as 
calculated for bicollinear (top) and collinear (bottom) Fe ar- 
rangements. Atom labels are as in Fig. [l] 



we have examined the iron or oxygen inserted systems (Fe 
or oxygen at A and B sites, as shown in Fig. [l] and Ta- 
bles U and [ill in addition to other modulations. Interest- 
ingly, the ground state magnetic structure of FeTe-Fe ap- 
pears to depend on the extra Fe content. As reported in 
Ref. [6], in Fei+^^Te, the favored magnetic structure is the 
bicollinear one as our calculations indicate for (5=0.125. 
For J = 0, the situation is reversed and the difference 
in total energy between the two phases with (5 = and 
5 = 0.125 is less than 0.2 eV (per supercell). The two 
magnetic states appear to be competing to be the ground 
state, with total energies that are quite close to one an- 
other. 

However, with oxygen insertion, this situation changes 
noticeably; with oxygen occupying either site A or site 
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FIG. 3: Comparisons of total energies of Fe-Te magnetic 
structures with various modulations. Note that the collinear 
structure is the more favorable one prior to going into a pos- 
sible superconducting state. 



B, the bicollinear arrangement becomes more stable with 
some rearrangements in the (magnitude and direction 
of) magnetic moments. However, with oxygen occupy- 
ing both A and B sites, the collinear structure (with re- 
duced/different moments) becomes the more stable one. 
We emphasize that this is an experimentally observed 
trend (Ref. [6 ). There is clearly no true antiferro- 
magnetism here; instead, we see some Fe moments re- 
duced by as much as 30-50% while others show smaller 
changes. We will label such magnetic states as damped 
SDW states. Our unit cell contains 8 Fe atoms and 
hence it is able to show the damping/oscillations bet- 
ter than what would be seen in a small (say Fe2Te2) 
cell; a larger supercell with more Fe atoms is likely to 
show the damping or oscillatory effects even better. In 
general, experiments show that static magnetism hinders 
superconductivity while spin fluctuations, which increase 
dramatically as the superconducting transition temper- 
ature Tc is approached, help or act as a precursor to it 
in the cuprates, pnictides and chalcogenides. Another 
point to note is the reduction in some Fe-site moments 
in the Feg-aTeg calculation with two vacancies in the 
collinear structure; however, in every case some notable 
local magnetic order is still present. 

We believe that, although the true, static antiferro- 
magnetism vanishes with the above modulations, local 
magnetic order plays a crucial role at least as a precur- 
sor to superconductivity. There is new evidence for such 
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TABLE I: Calculated collinear magnetic moments of Fe atoms 
for different configurations. Atom labels are as in Fig. [1] and 
initial magnetic moment direction changes are highlighted in 
bold face. 
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TABLE II: Calculated bicollinear magnetic moments of Fe 
atoms for different configurations. Atom labels are as in 
Fig.0 

behavior seen even in the cuprates. Recent RIXS (reso- 
nant inelastic X-ray scattering) experiments on the 123 
cuprate have revealed an intense peak around 1.7 eV en- 
ergy loss that is due to an optically forbidden d — d tran- 
sitions of unpaired holes from Cu2+ to other d-orbitals 
(Ref. [2D])- The significance of this work is that in over- 
doped, underdoped and superconducting samples, mag- 
netic excitations are seen around the same energy region. 

Total energies 

A comparison of variationally evaluated total ener- 
gies is one of the most useful and reliable outputs of 
a DFT-based calculation. Various modulation-induced 
structures were selected since, experimentally, these were 
identified as necessary precursors for FeTe-derived super- 



conductivity. Total energies shown in Fig. [3j for all the 
modulated structures show a very clear trend; i.e., the 
collinear-derived magnetic structure becomes the more 
favorable one prior to possible superconducting transi- 
tions. 

Here is a brief summary for the various structures: (a) 
Pure FegTeg: Collinear structure is more stable with Fe 
magnetic moments about 2iib per atom, (b) FegTeg: Bi- 
collinear is more stable with the extra Fe atom carrying a 
moment of about 2.5 fiB- (c) Oxygen interstitial: Insert- 
ing one oxygen atom either at the interstitial site A or B 
brings the total energies of the two magnetic structures 
close to one another and two inserted oxygens make it 
even closer. With three inserted oxygens, the collinear 
structure clearly becomes the more favorable, (d) Intro- 
duction of a single vacancy into the supercell induces a 
clear separation of the total energies with collinear one 
being more stable; nevertheless, a second vacancy brings 
their energies closer together, (e) Se doping makes the 
collinear structure more stable. 

While a one particle mean- field theory alone is unlikely 
to explain pairing related phenomena in moderately cor- 
related systems, the trends seen with oxygen insertions, 
Fe vacancies and Se doping appear to confirm several 
experimental findings, such as short-range magnetic or- 
der pi]. The true, static AFM state, with large Fe mo- 
ments, is lost due to oxygen insertions and other doping 
while some remnants of this lost state can be gleaned 
from the band calculations. However, without further 
assumptions and work, these results alone are unable to 
demonstrate a pairing scenario as to where the holes are 
created and their role due to these modulations. This is 
not surprising since it is well known that the local ap- 
proximations to exchange-correlation and mean-field ef- 
fects usually wipe them out. In addition, Fermi surface 
topology-based nesting mechanisms have become popu- 
lar and have been used to obtain various pairing mech- 
anisms. However, in view of the loss of long range an- 
tiferromagnetic order and the onset of short range in- 
commensurate order observed experimentally (Ref. [21]) 
and suggested by our supercell calculations, we present 
a different pairing scenario, based on a spatially local 
mechanism, in the following section. 

Charge and spin Pairing 

Here we introduce a many-body, Hubbard Hamiltonian 
where explicit fluctuations of paired holes can be system- 
atically studied in a grand canonical ensemble |22] |23] . 
Such calculations on a 8-site Betts cluster-based lattice 
(see Fig.|4| show that there are instabilities/fluctuations 
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FIG. 4: Negative charge gap and positive spin gap in the 8- 
site Betts-cluster-based lattice. The 8-site Betts cell, which 
is periodically repeated with edge vectors (2,2) and (2,-2), is 
also shown here. 



in charge and spin degrees of freedom under "suitable" 
conditions. These fluctuations exist between a back- 
ground AFM/SDW state and a two hole- or electron- 
doped state in a phase region identified as having a neg- 
ative charge gap favoring charge and spin pairing (at a 
critical value of doping). Although our recent work [24] 
refers to a single orbital model, we infer that it can be 
extended to two or multi-orbitals and still retain the neg- 
ative charge gap, as long as hopping strength among sim- 
ilar orbit als is greater than that for different ones. 

Depending on the value of the electron-electron 
Coulomb repulsion [/, there are certain regions in a 
3-dimensional T(temperature)-/i(chemical potential)-!! 
phase space where quantum critical points (QCPs) exist 
and drive these fluctuations. For the 8-site Betts cluster- 
based lattice, with periodic boundary conditions, such a 
QCP has been found dit U = Uc = 8.54 (in units of 
the hopping parameter in the Hubbard model), so that 
a charge instability exists in the region < U < Uc dit 
suitable (doping) ji = jic and T. This charge instabil- 
ity leads to fluctuations between a SDW/AFM state and 
a state with extra holes of paired (and unlike) spins at 
low temperature. Note that even for very small U (i.e., 
weakly correlated systems), a negative charge gap and a 
coherent (positive) spin gap can exist. Our grand conical 
ensemble-based studies show that these paired spins are 
stabilized by the (positive) spin gap up to a temperature 




FIG. 5: BicoUinear (Top) and collinear (bottom) magnetic 
structures in the plane containing Fe atoms. Arrows identify 
nearest and next nearest neighbors of the seventh atom (as 
marked) in the plane. 

Tg. Clearly, with extra holes, it is easier for the charge 
carriers to move between sites, as seen in the oxygen- 
doped case. 

An important point here is that such a scenario fa- 
vors a background of collinear-derived states discussed 
in the previous sections, rather than the bicollinear one. 
This is simply due to the differences in the neighboring 
spins which begin to appear at the next nearest neigh- 
bor (nnn) level of the two structures (see Fig. |5|. For an 
electron (or a hole) to hop through the lattice and form 
a pair, existence of (as many) unlike-spins as next near- 
est neighbors would be an asset (since like-spin neighbors 
would prevent hopping to that orbital) provided that U 
values are relatively small. This is what is seen in the 
collinear structure with 4 unlike-nnn-spins, compared to 
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2 unlike-nnn-spins in the bicollinear structure (while at 
the nearest neighbor level, there is no difference). Our 
calculated charge gap for the 8-site Betts-cluster-based 
lattice shows that the next nearest neighbor hopping, un- 
der certain conditions, can play a crucial, helpful role in 
the charge and spin pairing instability [24j . 



Summary 

Our band calculations, which are consistent with sev- 
eral experimental results, are used to identify states that 
act as precursors to superconductivity. In every case 
considered for the Fe-chalcogenide under consideration, 
a collinear-derived, damped SDW is predicted. A the- 
ory based on the Hubbard model, with a weak onsite 
Coulomb repulsion, is able to explain possible charge and 
spin fluctuations starting from these precursors. We be- 
lieve that this work lays the ground work for a common 
understanding of superconductivity in the chalcogenides, 
pnictides and cuprates. 
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